What do lattice baryonic susceptibilities tell us 
about quarks, diquarks and baryons at T > Tel 
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Lattice data on QCD thermodynamics, especially recent study of high order susceptibilities by UK- 
Bielefeld collaboration, have provided valuable information about matter properties around and 
above the critical temperature T c . In this work we tried to understand what physical picture would 
explain these numerical data. We found two scenarios which will do it: (i) a quark quasiparticle 
gas, with the effective mass which is strongly decreasing near the phase boundary into the QGP 
phase; or (ii) a picture including baryons at T > T c , with the mass rapidly increasing across the 
phase boundary toward QGP. We further provide several arguments in favor of the latter scenario, 
one of which is a natural continuity with the baryon gas picture at T < T c . 



I. INTRODUCTION 

QGP is experimentally studied via heavy ion collisions, 
at CERN SPS and last years at BNL RHIC collider, at 
temperatures reaching up to about T « 2T C . Success 
of hydrodynamical description [1] of observed collective 
flows have indicated, that all dissipative lengths are very 
short and thus the produced matter cannot be a weakly 
coupled gas but rather a near-perfect (small viscosity) liq- 
uid [2] . These features are further complemented by very 
high jet losses and robust heavy quark charm (equilibra- 
tion) observed, well beyond what pQCD predicted. As a 
result, a radically new picture of QGP at such tempera- 
tures is being developed, known as the strongly coupled 
quark-gluon plasma, or sQGP. 

It has been pointed out by Shuryak and Zahed [3] 
that the interaction seems to be strong enough to pre- 
serve the meson-like bound states above T c although in 
a strongly modified form. In particular, the lowest char- 
monium states J/ip,rj c are predicted to exist up to T as 
high as about (2.5 — 3)T C . These charmonium states were 
observed on the lattice [4] as peaks in spectral densities 
of the correlation functions, and they indeed seem to sur- 
vive till such high temperatures. 

It was further pointed out in the next paper by Shuryak 
and Zahed [5] that in the deconfined phase also multiple 
binary colored bound states should exist, in about the 
same T domain, since the interaction is about the same. 
To put the discussion below into proper perspective, they 
argued that there should be 3 categories of bound states, 
in decreasing robustness: (i) glueballs, (ii) ((73)3 and 
mesons qq; and (iii) (qg)e, diquarks and baryons. If the 
strength of the effective potential in qq states is counted 
as 1, the relative color Casimirs for categories (i) , (ii) and 
(iii) are 9/4, 9/8 ~ 1 and fts 1/2, respectively. In our 
recent work [6] we have extended the same approach to 
some many-body states. We found new 3-gluon configu- 
ration ggg belonging to category (i) , the polymeric chains 
q.g.g...gq of the category (ii) and diquarks and baryons 
in category (iii). 

The last two are the baryon number carrying states 



we will discuss in this work. Since these states belong to 
the third, most weakly bound category, they are natu- 
rally most vulnerable to uncertainties of the potential and 
their existence can be questioned. Besides, these states 
are relatively heavy: such states have not been included 
in [5] in pressure. 

The reason we will discuss them now is because they 
are more important at increasing baryonic chemical po- 
tential fi. Alternative way to look at the same thing 
is to consider higher derivatives over fj, at fi = 0: this 
way the role of such states is enhanced due to powers 
of their baryon number. At some point the diquarks 
and baryons should become noticeable in these quanti- 
ties even if their role in pressure is small: and this is 
precisely what we think happened in the lattice data of 
the UK-Biclcfeld collaboration (UKB) [7], especially in 
susceptibilities with 4 and 6 derivatives. 

In this work we will concentrate on the so called bary- 
onic susceptibilities part of the free energy, which can be 
singled out via derivatives over quark chemical potentials 
Mg = (M« + Mrf)/2 and fij = (/i u — /i<j)/2 calculated re- 
cently by the UKB. They use it in a context of Taylor 
expansion of the thermodynamical quantities in powers of 
baryonic chemical potential n/T* up to the order 0(/.i 6 ) 
of 2-flavor QCD, but we will not discuss this expansion 
per se and concentrate on (T-dependent) susceptibilities 
of the kind^: 



d n {T) 



d n (p/T 4 ) 



d(n/T)' 



i\c n {T) 



(1) 



fi=0 



for n = 2, 4, 6. (The odd ones vanish at /1 = by symme- 
try.) These data are shown in Fig.l and also below. The 



*We follow notations used in this work where ji is the chem- 
ical potential per quark, not per baryon. Thus the associated 
charge is B = 1 for a quark, B = 2 for a diquark and B = 3 
for a baryon. 

' Since we would not discuss any Taylor series in this work, 
we would prefer to leave out the factorials and thus discuss 
susceptibilities d n defined without them, not c n . 



1 



UKB also studied what they called isospin susceptibilities 
defined as 



dim 



d n {p/T 4 



= n!<£(T) (2) 



9WT)«- 2 a( A1/ /T)2 
and in a flavor diagonal-non-diagonal language there are 
<C - (d* + <)/4 , d u n d = (d n <)/4 . (3) 

Let us also mention another recent independent lattice 
studies on susceptibilities of 2-flavor QCD in [8], where 
they defined so-called nonlinear susceptibilities(NLS) 



Q(n u +n d )p 



(4) 



and evaluated their values on lattice. To see the connec- 
tion between these two approaches, we give the following 
relations 



71-2 



dT = T 



n—4 



2=0 



and 



jud 



(5) 



(6) 



While the two approaches are very closely related, their 
numerical results, however, are not quantitatively compa- 
rable, partly because they have used very different mass 
setup in the lattice calculation. We nevertheless empha- 
size that in a qualitative view both of them have found 
very similar and interesting patterns in those suscepti- 
bilities, especially for the 4th and the 6th, which are the 
central issues to be addressed in this paper. 

To set the stage, we start with the hadronic phase be- 
low T c . Here the relevant states are only the baryons 
with the baryon number (per quark) 3. Their spectrum 
is known at T — experimentally, and thus an obvi- 
ous question is: Can a simple resonance gas of known 
baryons explain the behavior of these susceptibilities be- 
low T c ? Indeed it is the CcLSG, clS shown by the dotted 
curves in Fig.l (obtained by including contributions of 
nucleon states from N(940) to N(1675) and A states from 
A(1232) to A(1700), for two-flavor theory one should not 
include strange baryons). No T— or ^— dependence of 
these masses is assumed, nor do we take into account 
the fact that lattice is dealing with non-massive quarks^: 
tuning these will shift the curves down a bit, making 
the agreement even better. So the susceptibilities in the 
hadronic phase, T < T c , can be described by the usual 
resonance gas of baryons. 



*In fact the input quark mass in these calculations is AT. 
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FIG. 1. The dotted lines correspond to a gas of baryonic 
resonances, the points with error bars are the susceptibilities 
d2(T), di(T), de{T) from lattice data (after removal of facto- 
rials in C2,4,e). 

The main issue to be discussed in this work is what 
these lattice data actually tell us about the nature of 
baryonic states above T c , and whether one can describe 
them with sQGP model or as well with some other model. 

Before we proceed to the argumentation in literature, 
let us remind standard thermodynamical expressions for 
massive fermions, which can be put in the following well 
known form 



JL = N M 



27T 2 T 2 



E 



(-) 



i+1 



I 2 



,lBn/T 



+e - lB »/ T )K 2 (lM/T) 



(7) 
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where B, M is the baryon number of the corresponding 
particle and its mass, N is the statistical weight and K 2 
is the Bessel function. § This form is very convenient for 
taking derivatives over /z, for example the first derivative, 
the baryon density is 

^=NB^f\tm (e iB, / T 

T3 ^2tt 2 T 2 ^ I [ 
1=1 l 

- e - lB ^ T )K 2 (lM/T) 

= NB Af[Bfj,/T, M/T] (8) 

where the function Af[x, y] is defined by these series. Note 
at this point we don't really consider mass as depending 
on jU so no extra derivatives against M appear. 

In a number of talks Karsch (and also a recent 
preprint) [9] have presented what we would refer to as 
a "naive" argument: the subsequent ratios 

dn+2/dn ~< B 2 > (9) 

are directly related to the squared baryon number of the 
constituents. The argument goes as follows: (a) For mas- 
sive particles with M >> T one can use the so called 
Boltzmann approximation, keeping only the first term in 
the sum above; (b) after that the fi— dependence fac- 
torizes, and thus each two derivatives over fi restore the 
same expression, modulo the factor B 2 . In the matter 
dominated by quark quasiparticles, or qg bound states, 
the r.h.s. would be 1, but it would instead be 4 or 9 for 
matter dominated by diquarks or baryons, respectively. 
The measured ratio d^jd 2 is <~ 10 at T < T c but at 
T > T c it rapidly drops and becomes close to 1. Com- 
paring it to the formula above Karsch concluded that at 
T > T c matter is a gas of some B = 1 objects, while the 
contribution of the B = 2 diquarks is strongly restricted. 

But if one looks closer at this argument, one finds it 
missing a lot of effects that should be there as well. For 
example, the next similar ratio d§ / d<± above T c is nowhere 
close to 1 but is in fact a large negative number <~ — 10 
which cannot be interpreted as a B 2 of anything. 

Furthermore, the idea that one can keep only the main 
term in the sum so that the /j, and T dependence can 
factorize, must be wrong by itself. The T-dependencc 
of d 2 (T), d±(T) and d§{T) is not at all similar: while 
d 2 (T) resembles the behavior of the pressure itself and 
can easily be interpreted as a transition from hadron to 
quark gas, the next one d\ (T) has a sharp maximum near 
T c , with even more complicated "wiggle" in the de(T). 



s If there are more than one species of particle we then sum 
over different species. Yet there will be particular concern 
when dealing with quasiparticles instead of particles where 
some background term may arise in the pressure, as will be 
discussed in later section. 



Another perspective on that issue can be made if one 
converts baryon number and isospin susceptibilities into 
flavor-diagonal (uu or dd) and flavor non-diagonal ud 
susceptibilities. The lattice data show that the second 
flavor- mixing derivatives are small** d 2 d /d 2 u << 1, but 
similar ratios for higher derivatives n=4,6 are not small 

<C7<C ~ i/2. 

Does it imply that the quark gas model is also inad- 
equate and should be excluded as well as the "bound 
state" gas? Or, if the argument is wrong, what exactly 
is missing? 

(i) Even if the Boltzmann approximation (keeping the 
first term in sum in (7)) may be good for pressure, it 
still fails for higher susceptibilities because the Z-th term 
has l[i in the exponent, and subsequent differentiation 
their role grows as l n . By the time one comes to the 
sixth derivative, these terms start canceling each other. 
In physics terms, this is a form of Fermi blocking effect 
not included in the simple Boltzmann approximation. 

(ii) The second to recognize is the fact that quasipar- 
ticles are not particles and their effective masses depend 
on matter parameters, such as T and especially fj,. Sub- 
sequent differentiation of this effective mass over /i would 
add powers of derivatives like 

d 2 M 

M" = —(T^ = 0) (10) 

to susceptibilities and to their ratios such as (9). Pro- 
vided those are large enough, they may completely inval- 
idate the naive interpretation of those ratios as baryon 
number squared. This was already pointed out by 
Bluhm, Kcmpfer and Soff [11], and we will refer to it 
below as the "BKS effect". The same is true for bound 
states such as baryons, and similar derivatives of their 
masses Mg(T) would play an important role below, 
(iii) the contribution of diquarks has been grossly overes- 
timated, while the contribution of baryons was not dis- 
cussed at all. We will show below that it may naturally 
explain the features seen in higher derivatives. 

The outline of the paper is as follows. In section II 
we will start with an "unconstrained" quark gas model, 
and will use the lattice data to extract the quasiparticle 
mass together with its dependence on matter, M(T,fi). 
We would not need to rely on perturbative arguments 
used by BKS [11] (since even their own fit leads to rather 
strong coupling at T w T c ). Furthermore, we will con- 
jecture possible relation between the T— and [i— depen- 
dences due to known shape of the phase boundary on the 
phase diagram. In section III we will further impose a 
number of constraints on quark mass, from other lattice 
data and also from confinement, a condition that there 
should not be any colored degrees of freedom at T < T c . 



"This is also the main point of the paper [10]. 
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We will conclude that these constraints basically make it 
impossible to ascribe the observed features of the data to 
the BKS effect. After that we will proceed to section IV 
in which we will discuss the contribution of diquark and 
baryons: here we will find good fits to the data satisfying 
all the needed constraints and nicely joining the baryon 
gas picture below T c . 



C?4 



d 3 (n B /T 3 ) 



dp, 3 
2tt 2 



/ 



dxx 



n 



^ (3) (e ) 



+3n 2 F^(e ) 



mo 



,d 2 m 



dfl 



7,2 \ ll =o> 



(13) 



II. MODEL I: A QUARK GAS WITH AN 
UNCONSTRAINED MASS M(T, fi) 

The idea to use thermodynamical quantities calculated 
on the lattice to fit the mass parameters of quasiparticles 
is by itself quite old. For example, Levai and Heinz [12] 
have used the data on p(T) for determination of quark 
and gluon effective masses M(T)^ . 

One well known problem with quasiparticle gas models 
is that the derivatives over T and fi upset thermodynam- 
ical consistency between gas-like expressions for different 
thermodynamical quantities. Only one of them can be 
assumed to have a simple additive form over quasiparti- 
cles: then there is no freedom left and all other quantities 
can be calculated from it by thermodynamics. Thus only 
one "primary" expression can be additive, while others 
will have extra "derivative" terms complementing simple 
gas formulae. 

Following conventions of the BKS paper, we will use 
as such "primary" expression that for the baryon num- 
ber density (8). The expressions for pressure and energy 
density would then be corrected by some T, /i-dependent 
"bag terms". Higher derivatives terms d n will be cal- 
culated by differentiating (8) n — 1 times. To be more 
specific, we explicitly give the baryon number density for 
this quark gas model 



n B 
T 3 



d(p/T 4 ) 
d(fx/T) 



-^t- I dxx 2 n [F(e — njl) — F(e + njl)} 



(11) 



Here g = N s * N c * Nf = 12 is the degeneracy factor for 
quarks in the two-flavor case, n is the baryon quantum 
number of quark which is defined here to be n = 1 by 
setting /i to be the quark chemical potential, fi = fi/T is 
made to be dimensionless, and e = V x 2 + fa with m = 
M/T. And finally we have introduced Fermi distribution 
function F(y) — Starting from (11) the explicit 

formulae for d 2 , d±, are given to be: 



d 2 



d(n B /T 3 ) 



djl 



-|L J dxx 2 n 2 F^(e ) (12) 



"it was not as direct as our approach below, because one 
cannot get 2 functions out of one without assumptions. 



d 6 = 



d 5 (n B /T 3 ) 



dfl 
2tt 2 



r,5 



A»=0 



dxx 



too ,d 2 fn I 



+ 10„-F«>(«,)^( 



eo K d[i 2 
to 2 ,d 2 fa 



-J 



+ 



3x 2 



x 2 + toq 2 dfi 



,0 'to, >: 



(14) 



In the above equations we have used e = \/ x 2 + m 2 
and too (T) = M(T,fj, = 0)/T, and also F^(y) means 
the zth derivative of the function F(y). 

The model used in the BKS paper assumes some Hard 
Thermal Loop based perturbative form for the T, /i- 
dependent mass with the coupling g 2 (T, /i) running in 
a complicated fashion fitted to reproduce the susceptibil- 
ities we discuss in this work. However, we do not see why 
any assumptions about the mass dependence are actually 
needed^ at this point. 

We thus suggest a generalization of what was done in 
[11]. Assuming a simple ideal gas model of quark quasi- 
particles, one has their mass to be the only input needed. 
With the lattice data on ^(T 1 ), d±(T) and de(T) used as 
input, one can simply solve for the three functions of T 
which would ideally fit them: we have chosen those to 
be: (i) the quark mass M(T,/j, = 0); and its two lowest 
non-zero §§ derivatives over \i (ii) M" = $pf(T,n = 0) 

and (iii) M"" = 3jjr(T,/i = 0). With these at hand, of 
course, we are able to develop the Taylor's expansion for 
quark mass as a function of | ^ | < 1 : 



M(T,£) = M(T,/i = 0)+i^ 



"There is of course no reason to trust any perturbative for- 
mula near T c at all, where the coupling becomes as strong as 
it was found by BKS themselves. 

^The quasiparticle masses and other quantities obviously 
can depends only quadratically on fx because of ft — > —fi sym- 
metry based on CP invariance. 
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i a 4 M, 



(15) 



unknown term is 



d m 



which could be solved out 
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FIG. 2. Quark quasiparticle mass and its second and fourth 
derivatives over \i as a function of temperature T, extracted 
from lattice data for susceptibilities. There are two sets of 
points in each figure that are obtained from C2, C4, C6 and from 
c 2i c 4) c 6 respectively. In the top figure for quark mass, we also 
plotted the two points with error bars measured by lattice via 
propagator, and the mass given by (20) as well, (the dashed 
line) . 

The procedure is iterative: First we used c^T) data to 
solve for the mass to as unknown. Then we go to a, the 
equation of which includes both to and §J?|^ =0 , but 
since we have already solved too from c 2 now the only 



from lattice results of C4. 
from cq with too and 



d£ 4 l M =0 

already being solved from 



finally we can obtain 

C2 and C4. The results for these three functions are shown 
in Fig. 2. (The error bars in too are determined from 

r) 2 ~ I 

uncertainty in C2- While for ^j?|^_ the errors should 
come from both C4 and too, the error bars in the figure 
only include those from C4, and also for ^pr| M _ the error 
bars solely include that originated from c§.) 

As an independent check, we have also extracted the 
same three quantities, M(T,fi = 0), ^0-(T,n = 0) and 
^-{T,n = 0) from the lattice data set for c T 2 (T), c{(T) 

and Cg(T) from [7] by the same strategy (but starting 
with isospin densities). 

The results are shown in Fig. 2. As can be seen, two sets 
of parameters we extracted from both data sets are well 
consistent with each other at T > T c , while for T < T c 
they do not agree. It is a good feature, as the quark gas 
model is not supposed to work there, in the domain of 
the baryon resonance gas. 

Let us summarize these results. The most important 
lessons are: (i) the mass M(T) strongly increases when 
cooling down toward the critical point T c ; (ii) Large and 
negative ^M(T, fi = 0) close to T c ; (iii) The 4-th deriva- 
tive is positive: so this decrease of the mass due to the 
2nd derivative will stop at about u/T ~ 1, see (15). 

The first two points are the trends already emphasized 
by BKS [11]. In their approach these two features are 
related with each other because of the assumed pertur- 
bative origin of the effective quark mass: 



M 



9 2 (T, V)T 2 (l + N f/G + ^ E tij ( 16 ) 



where the sum runs over all flavors /. Ignoring for a 
moment a (rather complicated) running of the coupling, 
the BKS mass is thus constant at the particular ellipsoids 
in the T — n plane, thus the derivatives over T and u arc 
related. 

We would like to propose another reasoning that leads 
to similar effect, but is free from perturbative assump- 
tions. Its idea can be described as follows: the quark 
mass should be getting large not only near the critical 
point T — > T c , u = 0, but near the whole critical line, 
at all a. It is needed to ensure that quark degrees of 
freedom do not contribute in the confined phase, at any 
fi. 

The critical line at nonzero fj, is schematically shown 
in Fig. 3, its shape at not-too-large \xjT can be described 
by an ellipsoid, or an unit circle, if the units are chosen 
appropriately. One may further think that the mass de- 
pendence on the radial coordinate R on such a plot is 
much more important than on the angular one <j> since 
the "lines of constant mass" should be nearly parallel to 
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the critical line, at least in its vicinity where the discussed 
effect takes place. 



T/Tc 




A, 4> 



|Ll/|Ll r 



1 r~r~ c 

FIG. 3. In the plane of temperature T-baryonic chemical 
potential /u, both appropriately normalized, the phase bound- 
ary looks like a part of a circle. (At least for the part marked 
by the solid line, studied well at SPS and RHIC, with quite 
well established chemical freezout. The dashed line is a con- 
tinuation of the freezout line where its association with the 
critical line is questionable.) The polar coordinates to be used 
are the radial distance R and the angle 4>. 

So, the proposed extension of the T-dependence of the 
mass to its /x-dependence is based on a substitution 



M(T,ii = 0)^M(R(T,n)) 



R 2 = t^ + ^ 



(17) 



(18) 



We have introduced here a new parameter fj, c : its value 
can be readily obtained from the experimental freezeout 
curve measured in heavy ion collisions at small /i, be- 
lieved to represent the critical line. If so, the value of 
this parameter is 



[i c /T c = 1.7 



(19) 



which is quite different from the value given by "pertur- 

bative scaling" (16): \J (I + ^-)k 2 w 3.63 which is not 
supposed to work in the non-perturbative regime near 



contradicts direct lattice measurements of the quasipar- 
ticle masses ; (ii) It implies that quark degrees of free- 
dom still significantly contribute in the confining phase 
at T < T c . In this section we will show what happens if 
one tries to modify the unconstrained model to make it 
compatible with both. 

One feature of the Model I is the relatively light quark 
mass M(T, \i = 0) in region 1 — 2T C ranging from about 
1.7T C to 2.2T C . Such mass conflicts with another lattice 
data about quark quasiparticle mass at 1.5T C and 3T C , see 
[13] which are m q /T = 3.9 ± 0.2 and m q /T = 1.7 ± 0.1, 
respectively, and are shown in Fig. 2 by two crosses with 
the error bars. Although these results are based on only 
one paper and have not been systematically studied by 
other lattice groups so far, they nevertheless represent di- 
rect measurements from the quark propagators. Further- 
more, such large masses correspond to the inter-particle 
potentials at large distances measured in separate lattice 
study [14]. 

Although the mass extracted via the Model I grows 
toward T c , this effect is still not robust enough to make 
quark contribution near-zero (or negligible) at T — T c . 
(In fact, BKS proceeded to fit equally well some region 
below T c .) This is unacceptable, since we know that there 
are no propagating quark degrees of freedom in the con- 
fining phase. 

Both these reasons force us to reconsider Model I, ba- 
sically by increasing the quark mass significantly to meet 
both constraints. This can be achieved by a quark mass 
formula similar to that used in [5] 



M(T) 



0.9 

T — 1 



+ 3.45 + 0.4T 



(20) 



with all units in proper powers of T c (This and subse- 
quent mass formula would then be generalized to finite 
/i according to (17) . The coefficients are chosen so that 
the curve goes through the two lattice-measured points 
for quark mass at T = 1.5, 3.0 T c , see the dashed line in 
Fig.2. 

We show what happens then to the susceptibilities, see 
the medium-thickness solid lines in Fig. 7. In short, good 
description of c 2 (T) is definitely ruined The issue is the 
same as for pressure in [5] , and perhaps can be cured by 
qg and other bound states. But this is not the only prob- 
lem of the constrained model: although it can produce a 
peak in di{T) and a "wiggle" in d§{T), given large enough 
derivatives over /i, those get displaced toward larger T as 
compared to the data. It is an inevitable consequence 
of the second constraint, insisting that quark effect be 
effectively zero at T c .*** 



III. MODEL II: THE CONSTRAINED QUARK 
GAS 

The "unconstrained Model I" discussed above, al- 
though consistent with both data sets d n {T) , d^(T) , is 
unfortunately unacceptable, for two main reasons: (i) It 



***The very heavy mass due to the constraints significantly 
decrease quark contribution to thermodynamics and hence 
disfavor quark-only model, yet on the other hand, it strongly 
favors the formation of bound states. 
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Completing our discussion of purely quark models, we 
now proceed to the possible role of their bound states, 
diquarks and baryons. 

IV. THE EFFECT OF DIQUARKS AND 
BARYONS 

We will now proceed to contributions of the bound 
states to the baryonic susceptibilities. Let us remind the 
reader that the particular reason to focus on diquarks 
and especially baryons is that the role of diquarks and 
baryons relative to quarks grows with fi because of their 
larger baryon charges. Alternatively, their contribution 
to the susceptibility d n grows exponentially with n: by 
the factors 2™ for the diquarks and and 3™ for baryons. 
For example, the contribution of N, A is enhanced by a 
factor 81 for <i 4 and 729 for d 6 relative to pressure: with 
estimates of pressure given above one may then expect to 
see their contribution there. On the other hand, for lower 
derivative (fe we expect quark-gluon bound states, which 
are more numerous and more tightly bound, to contribute 
significantly. We summarized all bound states, together 
with their multiplicities, in the Tabled. 

(In passing, let us comment about the numbers in the 
real world with strangeness, Nf = 3. The number of 
diquark flavor states is increased to be 3 times larger, for 
baryons the total spin-flavor multiplicity increases from 
4+16=20 to 56(an octet J = 1/2 and a decuplet J = 3/2) 
which is roughly enhanced by 3 times, so the numbers 
both diquarks and baryons states are increased by the 
factor 3. The quark number increases as 3/2, so the 
overall enhancement of the ratios we will discuss below 
from N f = 2 to N f = 3 is the factor 2.) 

TABLE I. Summary of states with baryon number at 
T > T c studied in this paper. 



state spin flavor color multiplicity 

q 2 2 3 12 

(93)3 4 2 3 24 

(qg)e 4 2 6 48 

(qq)t I= ° 113 3 

(qq)i- 1 - 1 3 3 3 27 

N 2 2 1 4 

A 4 4 1 16 



Quark-gluon bound states: Before we proceed to 
actual calculation, let us make simple estimates of the rel- 
ative weighing in pressure of quark-gluon bound states. 
The 2-body states qg are thermodynamically suppressed 
by additional Boltzmann factor, exp(—M/T) ~ 0.02 — 
0.04 (by including their considerable binding). However, 
due to their relatively large multiplicity (6 times the num- 
ber of the quark states) they contribute to the pressure 
and susceptibilities at the level of about 1/10 or more. 

To get more quantitative answer one has to know the 
binding energy of these states. While the binding of the 
category three states (qg)6 can be reasonably neglected, 
the category two (93)3 states have considerable binding 
at the same order as meson states. The potential model 
calculations in [6] lead to (qg)3 binding up to \5E\/T rts 
1.4 at T = T c , which means their contribution increases 
relative to simple estimate above by extra factor 2-3. 

However there are many reason to doubt that close 
to T c this calculation can be trusted quantitatively. In 
particular, the potential used is measured on the lattice 
for static charges only, and the corresponding calcula- 
tions are supposed to be reliable only when the binding 
is small: near T c more complicated dynamics beyond the 
potential model will contribute as well. 




1 1.5 2 2.5 3 

T/T c 

FIG. 4. Comparison of susceptibilities from quark-gluon 
states with two limiting case, zero-binding ((^2,4,6) and 
"full-compensation" binding (£>2,4,6)- 

Let us thus just suggest an upper limit for the qg 
states' contribution. Since the qg states are colored, 
they should gets infinitely heavy at T c , together with all 
other colored states. Furthermore, as (the more tightly 
bound) qgs states have the total charge of one quark, 
their mass should not be smaller than that of one charge 
\SE\ < M q (T w T c ). So we expect M(qg) to interpolate 
between M q + M g w 2M q at zero binding to a a single M q 
at T w T c ). The contribution of these states to suscepti- 
bilities in the two limiting cases, namely the zero-binding 
case (labeled in figure by di) and "full-compensation" 
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binding case (labeled by Di), are shown in Fig. 4. We 
conclude that large uncertainty, of the order of factor 3, 
remains in the contributions of such states. 

(These results are calculated with (qg)e always having 
twice quark mass and melting at 1.4T C while with (qg)s 
having twice quark mass in the former case and the same 
mass as quark in the latter, both melting at 2.1T C . The 
actual contribution of quark-gluon bound states should 
be somewhere in between, near to Di around T c while 
rapidly decreasing to di for higher temperature.) 

One may also ask what is the contribution of the var- 
ious polymer-like qg states qg, qgg, qggg, • • • which, ac- 
cording to [6], has the same binding energy per bond. 
The effect of these states can be easily evaluated via a 
geometric series: the resulting enhancement factor is 



f, 



1 



polymers 



1 - 6exp[(\SE\ - M g )/T] 



(21) 



where 6 is the color and spin degeneracy added by each 
link. For small binding this is just a few percent cor- 
rection, but if it may get to be strong enough to drive 
the denominator toward zero, a total "polymerization" 
of sQGP would occur. 

Diquarks: for Nf — 2 gauge theory corresponding to 
the UKB data at hand there is only one attractive di- 
quark channel, the antisymmetric color triplet*** (33)3- 
Because of Fermi statistics, it means that the product of 
spin and flavor should be symmetric, and thus there are 
two options: (i) spin-0 isospin-0 ud diquark (qq)^ 1 ^ , 
and (ii) spin-1 isospin-1 one (qq)^ =I=1 . These are the di- 
quarks which are familiar in hadronic spectroscopy, the 
former appears inside the N, the latter inside A (octet 
and decuplet members, for 3 flavors). The lesson from 
this spectroscopy (at T=0, of course) is that while the 
former is well bound, by about 300 MeV, the latter is 
not. In view of the rather marginal character of diquark 
binding, we expect only the former one able to be se- 
riously considered as bound state above T c . Neverthe- 
less to confirm the point that diquarks will not play any 
role in all susceptibilities measured, we include both of 
them in calculation of Fig. 7. If we only use antisymmet- 
ric states, then the contribution will be reduced to only 
1/10 of that. 

The diquark-to-quark pressure ratio can be estimated 
as following: 



(99)3 



12 



27 e Xpi "- M+m )^ 



(22) 



where the binding \SE\ is negligible (actually only the 
3 (qq)i =I=0 states are very likely bound)***. At small 



/i where the data under consideration are calculated, 
M/T w 5 and their contribution is at few percent level, 
negligible compared to uncertainties. 

Baryons: as we found in [6] they are bound till about 
T = 1.6T C . In the 2-flavor theory they are the N, A 3- 
quark states. Only the s-wave basic states survive above 
T c , while all other resonances (used in the first section at 
T < T c ) which are orbital or radial excitations of N, A 
families are "melted" . 

The baryons are also numerous (20) but the suppres- 
sion factor due to mass is much smaller 



(ggg) „ 20 2 M -2M+|^| 3/2 
q ~ 12 { T ' 



(23) 



Near the "endpoint" of baryons with zero binding (which 
according to [6] is at T = l.QTc) their mass is 3M q , 
expected to be in the range of 2.5- 3 GeV. As it is an order 
of magnitude larger than T, one would not expected to 
contribute to pressure etc. 

However, unlike the quark, quark-gluon and diquarks 
(which after all are colored objects existing only above 
T c ), N, A baryons are colorless and thus survive on both 
sides of the boundary of (a continuous) phase transi- 
tion (a crossover, more accurately), thus the masses of 
baryons at T — > Tc are expected to join continuously to 
their known values at lower T. This of course implies 
that the binding energy near T c gets very large due to 
some deeper yet poorly-known mechanism, and the po- 
tential model used in [6] to evaluate this binding will not 
be applicable. The situation is basically the same as with 
mesons: as emphasized in [5] the pion mass must (by def- 
inition of chiral breaking) vanish (in the chiral limit) at 
T — > Tc, which potential model also cannot reproduce, 
we will use below the following parameterization (in T c 
units) 

M N = 9.5 + 4.6 * tanh[3.8 * (T — 1.4)] (24) 



M A = 10.25 + 3.85 * tanh[3.8 * (T - 1.4)] (25) 

interpolating between the nucleon and A vacuum masses 
at low T, while approaching the same value 3M q at high 
temperatures. We plot it in Fig. 5, together with the 
masses of various other states to be used in later in Fig. 7 
for susceptibilities. The main feature is fundamentally 
again enforcing confinement: when going from QGP side 
toward T c , all colored degrees of freedom get extremely 
heavy and drop out from system, while all colorless de- 
grees of freedom get more tightly bound and eventually 
dominate. 



***In the talks [9] Karsch mentioned "hundreds of di- 
quark bound states" in our model: it must be some 
misunderstanding. 

^The last factor comes from M 3//2 in the pre-exponent, 
originated from momentum integral. 
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0.5 1 1.5 2 2.5 3 

T/T c 



0.8 1 1.2 1.4 1.6 1.8 2 

T/T c 



FIG. 5. Masses of various states studied in this work. The 
thin solid line is for quark and the dashed line is twice quark 
mass which is roughly for quark-gluon and diquark. The lower 
thick solid line is for nucleon states and the upper one for A 
states. These masses are used for calculation of Fig.7. 

We now approach the central point of the paper: the 
baryon contribution provides a natural interpretation 
of the structures observed in susceptibilities measured 
on the lattice, the large peak near T c in d±{T) and a 
more complicated "wiggle" structure is seen in d$(T). 
This happens because the expected mass dependence of 
baryons on T, fj,, shown in Fig. 5, should have a character- 
istic shape with an inflection point, separating the region 
in which the second derivative Mg is negative (above T c ) 
and positive (below T c ). That is why the contributions 
of the baryons to d 6 show a "wiggle" as seen from the 
corresponding curves in Fig.7. Note also, that there is a 
less pronounced wiggle of the same origin in baryonic d^: 
we think its negative part is the reason why the qg and 
qq contributions above T c can get compensated and by 
coincidence the d^jdi ratio gets close to 1 there. 

One additional argument for baryonic nature of the 
structures seen in di,d§ is the following one. Each 
derivative over n q leads to factor 3, so 2 of them give 
9. If instead one has two derivatives over \ii the fac- 
tor obtained is (2/ 3 ) 2 , which is 1 for p, n, A+,A° and 
9 for A ++ ,A~. As a result, if one ignores the mass 
difference between these states, one finds that baryonic 
contribution to both should have the ratio d^/dn = 
(1/9) * (4/20) + (1 + 1/9) (8/20) = .467, while this ratio 
should be 1 for ideal quark gas. The actual ratio of these 
quantities according to UKB data are shown in Fig. 6. 
We see near T c the data obviously favor the existence of 
baryons, especially for dl/d & , and the quark asymptotic 
end is arrived at about 1.4T C for d\jd^ while only after 
1.8T C for d^/de. These evidences strengthen the necessity 
of baryonic interpretation of the higher susceptibilities. 



FIG. 6. The susceptibilities ratios d\jd^ (the thin solid) 
and dg/de, (the thick solid). The dashed lines correspond to 
ideal quark gas (upper) and ideal baryonic gas (lower). 

Taking everything together, including quarks, 
quark-gluons, diquarks and baryons, we arrived at sum- 
mary plots shown in Fig.7. We repeat that all masses 
used are as shown in Fig. 5 and their \i— dependence is 
introduced in the same way according to (17). 

( The bound states' "endpoints" are set to be 2.1T C for 
{qg)z quark-gluons, 1.4T C for (qg)e, 1.4T C for diquarks, 
and 1.6T C for baryons, according to [6]. The gradual 
removal near melting point is done by similar means as 
in [5]. The results are shown in Fig.7, where the overall 
values as well as the contributions of each kind of states 
are all present.) 

Let's focus on the T > T c side. The conclusions are: 

(i) as expected the diquark contribution is negligible 
for all three quantities even after including the suspect 
(qq)i =I=1 states, but it is clearly growing as getting to 
higher derivatives; 

(ii) For d 2 quark provides main contribution, and we em- 
phasize the fitting will be much better if we include the 
large binding of qg states near T c . We have shown above 
that large uncertainty in its binding, including polymers, 
would allow for a good fit here, which we decided not to 
do. 

(iii) In g?4 it is precisely the baryons that produce the 
desired large peak near T c till about 1.3T C where quarks 
become important; 

(iv) The baryons' contribution extremely dominant the 
behavior of dg, especially the "wiggle" shape. 

We conclude that two prominent structures, a peak in 
cLi(T) and a "wiggle" in d§(T) are naturally reproduced 
by baryons. 
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FIG. 7. The contributions of different states to (a) c?2, (b) 
di and (c) ^6, as well as the summed total values. The thick- 
est solid lines are for taking all together, while the medium 
solid lines for quark, the thin solid lines for baryon, the dot- 
ted lines for quark-gluon, and the dashed lines for diquark, 
respectively. 



V. SUMMARY 

In one sentence, the main lessons from the UKB sus- 
ceptibilities is that the baryons NA do survive the QCD 
phase transition, but are rapidly becoming quite heavy 
across it. 

More generally, the discussed data set on the baryonic 
and isospin susceptibilities at T > T c can be described 
in two different scenarios, (i) The first is a quark quasi- 
particle gas, with the effective mass which is strongly 
decreasing near the phase boundary into the QGP phase; 
(ii) the second is a picture including baryons with the 
mass rapidly increasing across the phase boundary to- 
ward QGP, to about 3M q . 

The first scenario was already pointed out by BKS [11], 
while our discussion makes it a bit more general. Its at- 
tractive features notwithstanding, it suggests the values 
of the mass not large enough to accommodate the exist- 
ing constraints from other lattice measurements. We also 
think it is not possible to have quark degrees of freedom 
in hadronic phase. Thus we conclude that success of such 
scenario is unlikely. 

The second scenario, based on baryons, can provide 
another explanation of the main features of the data, 
namely the observed peak in d±(T) and a "wiggle" in 
de(T). It also naturally explains the flavor-changing 
d% d ,dQ d 7 which are not small relative to flavor-diagonal 
ones. Last but not least, this scenario provides a desired 
continuity to the baryon resonance gas picture at T < T c . 

Although the susceptibilities d n (T) we used in this 
work are highly sensitive tools, they are quite indirect. 
Thermodynamical observable in general cannot tell the 
difference between "melting" baryons (getting unbound) 
and baryons remaining well bound but just getting too 
heavy: in both cases all one finds is that their contribu- 
tion to thermodynamics effectively disappears. Besides, 
the ideal gas models used in these studies are probably 
too naive to claim really quantitative description of the 
data. One should instead study directly the spectral den- 
sities of the correlators of the appropriate baryonic cur- 
rents (qqq) and see if there are baryonic peaks there, like 
what has been done for charmonium and light mesonic 
channels. Only such direct measurements would tell us 
which scenario is the correct one. 

Speaking about experimental confirmation of the 
"bound state" scenario, we think the best chance could 
be observation of the vector mesons. As described in de- 
tail in [15], vector mesons p 1 lu, <fr are expected to become 
heavy near their disappearence point, like the baryons 
discussed above, reaching the mass w 2M q — 1.5— 2 GeV. 
The next generation of RHIC dilepton experiments have 
a chance to see if this is indeed what is happening in 
QGP. 
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